Abstract 



We consider small solutions of some vibrating mechanical systems with 
smooth non-linearities for which we provide an approximate solution by 
using double scale technique; a rigorous proof of convergence of the double 
scale method is included; for the forced response, a stability result is 
needed in order to prove convergence in a neighbourhood of a primary 
resonance. 
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1 Introduction 

In this work we look for an asymptotic expansion of small periodic solutions 
of free vibrations of a discrete structure without damping and with local non 
linearity; then the same system with light damping and a periodic forcing with 
frequency close to a frequency of the free system is considered (primary reso- 
nance). For a small solution, we recover a behavior with some similarity with 
the linear case; in particular the amplitude of the forced response reaches a 
local maximum at the frequency of the free response. On the other hand the 
frequency of the free response is amplitude dependent and the superposition 
principle does not apply. The work of Lyapunov [1] is often cited as a basis for 
the existence of periodic solutions which tends towards linear normal modes as 
amplitudes tend to zero; the proof of this paper uses the hypothesis of analycity 
of the non linearity involved in the differential system. In [2], we addressed 
the case of a non linearity which is only lipschitzian and we prove existence of 
periodic solutions with a constructive proof; in this case the result of Lyapunov 
obviously may not be applied. Non-linearity of oscillations is a classical theme in 
theoretical physics, for example at master level, see [3] in Russian or its English 
or French translation in [4j [5] . 

Asymptotic expansions have been used for a long time; such methods are 
introduced in the famous memoir of Poincare [5] ; a general book on asymptotic 
methods is [7j with french and English translations [8j [9] ; introductory material 
is in [10] , [TT] ; a detailed account of the averaging method with precise proofs 
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of convergence may be found in [12 ; an analysis of several methods including 
multiple scale expansion may be found in |13j ; the case of vibrations with uni- 
lateral springs have been presented in [H OS [IB] , [TH ITS] |T§1 12T H |2"T] : in [22 ] a 
numerical approach for large solutions of piecewise linear systems is proposed. 
A review paper for so called "non linear normal modes" may be found in [23] : 
it includes numerous papers published by the mechanical community; several 
application fields have been addressed by the mechanical community; for exam- 
ple in |24j "nonlinear vibro-absorption problem, the cylindrical shell nonlinear 
dynamics and the vehicle suspension nonlinear dynamics are analyzed" . 

In the mechanical engineering community the validity of the expansions is 
assumed to hold; however, this is not straightforward as this kind of expansion 
is not a standard series expansion and the expansion is usually not valid for all 
time; for example, this point has been raised in |25j . If the averaging method 
was carefully analyzed as indicated above, it seems not to be the case for the 
multiple scale method, the expansion of which is often compared to the one 
obtained by the averaging method. 

Here in a first stage we consider small solutions of a system with smooth 
non-linearities for which we provide an approximate solution by using double 
scale technique; a rigorous proof of convergence of the double scale method is 
included; for the forced response, a stability result is needed in order to prove 
convergence. As an introduction, the next section addresses the one degree 
of freedom case while the following one considers many degrees of freedom; 
for free vibrations we find solutions close to a linear normal mode (so called 
non linear normal modes) and for forced vibrations, we describe the response 
for forcing frequency close to a free vibration frequency. Preliminary versions 
of these results may be found in [26] and have been presented in conferences 
[2"T1 |2"5] ; related results have been presented in [33] . Triple scale expansions is 
in preparation by N. Ben Brahim 30 . In a forthcoming paper, the non-smooth 
case will be considered as well as a numerical algorithm based on the fixed point 
method used in [2]. Such vibrating systems linked to a bar generate acoustic 
waves; this point will be studied in an other work. 

2 One degree of freedom, strong cubic non lin- 
earity 

In this section, we consider the case of a mass attached to a spring; in the case 
of a stress-strain law of the form n = ku + mcu 2 + mu 3 , we find no shift of 
frequency at first order, so here we concentrate on a stress-strain law with a 
stronger cubic non linearity: 

n = ku + mcu 2 + m— it 3 , 
e 

where e is a small parameter which is also involved in the size of the solution as 
in previous paragraph; the choice of this scaling provides frequencies which are 
amplitude dependent. 
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2.1 Free vibration, double scale expansion up to first or- 
der 

Using second Newton law, free vibrations of a mass attached to such a spring 
are governed by: 

u + uj 2 u + cu 2 + — = 0, (1) 
We look for a small solution with a double scale for time; we set 

T = cut, Ti = et, (2) 
so with Doit = -Diit = we obtain 

dlt „ „ d 2 7i , o „ „ n n 

■ = u;L> w + erJi-u, —— = u 2 Diu + 2euiD D 1 u + e z D{u (3) 
at at 2, 

and we look for a small solution with initial data 

u(0) = ea + o(e) and u(0) = o(e); we use the ansatz 



so we have: 



and 

d 2 u 
dt 2 

with 



u = eu 1 {T ,T 1 ) + e z r{T Q ,T 1 ,e), (4) 

^ = e[wD ui + e-Dxtti] + e 2 [u}D r + eDir] (5) 

ew 2 L> 2 ui + e 2 [2o;£>o-DiWi + ^ 2 £>;H + £ 3 [£>?«i + ZV], (6) 

1 /d 2 r \ 

V * r = Z ( ^2 ~ w D o r ) = 2uD D ir + eDfr. (7) 



c 



We plug expansions (Q}, © into ([T]); by identifying the powers of e in the ex- 
pansion of equation ([l} , we obtain: 



uj 2 {D 2 Ui +m) = 0, 
(D 2 r + r) = %, with 



(8) 



S2 — — 7 



c(eui + e 2 r) 2 + -(eu\ + e 2 rf 
e 



2LjD DiUi-eR,(ui,r,e), (9) 



with 

K = D\ui + V 2 r; (10) 

we can manipulate to obtain: 

S 2 = - [cu\ + du\ + 2ujD q D iUi + eR( Ul ,r, e)] , (11) 

where 

R(ui,r,e)= \jl + 2cuir + 3du\r + ep(ux,r, e)] , (12) 
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with a polynomial p{u\, r, e) = cr 2 + idu\r 2 + edr . 

We set 0(T o ,Ti) = T + /3(Ti) noticing £> o = 1, Drf = D x j3; we solve 
equation flSJ) with: 

ui = o(Tx) cos(0) (13) 

and we obtain 

— rn 2 dn 3 

S 2 = -^-(1 + cos(20)) - — (cos(30) + 3cos(0)) + 

2oj(D ia sm{8) + aD 1 pcos(9j) - cR(ui, r, e); (14) 
we gather terms at angular frequency 1: 
da 3 

S 2 = — — 3cos(0) + 2w[L> 1 asin(0) + aDifi cos{9)] + - eR(y,i,r,e), (15) 
where 

Si = ^(1 + cos(20)) - ^ cos(30). (16) 

By imposing 

da 3 

£>ia = and 2ujaDij3 = 3^—, so that 

da 2 

a = a Q , = j3 Ti with O = 3— Ti, (17) 

we have canceled the constant involved in /3o with the choice of zero initial 
condition and we get that S 2 = S% — eR(ui, r, e) no longer contains any term at 
frequency 1. 

In order to show that r is bounded, after eliminating terms at angular fre- 
quency 1, we go back to the t variable in the second equation ([8]). 

r + uj 2 r=-^, with (18) 

id 

S 2 = Sl (t,e)-eR(u 1 ,r,e), where (19) 

dcji 

S*(t, e) = — ^— [1 + cos(2(u;i + j3(et)))] - — cos(3(o;i + f3(et))) (20) 

— co? d(X^ 

-(1 + cos(2(u;i + /3 a etj) — (cos(3(u;i + /3 et))) , (21) 



2 v ™ /, 4 

with ^(ui, r, e) = R(ui,r, e) - V 2 r, (22) 

in which the remainder R is expressed with variable t. 

Proposition 2.1. There exists 7 > swc/i £/ia£ /or all t < t e = - , the solution 
°f (EPt with u(0) = eao + o(e), w(0) = o(e), satisfies the following expansion 

u(t) = eao cos(i/ e t) + e 2 r(e,i), 
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where 

„ da 2 ,„„. 
^ £ = uj + 3e— (23) 

and r is uniformly bounded in C 2 (0,i c ) . 

Proof. Let us use lemma HOI with equation (|18p ; set 5 = S^; as we have enforced 
(fTTj) . it is a periodic bounded function orthogonal to e ±!t , it satisfies lemma 
hypothesis; similarly set g = R; it is a polynomial in variable r with coefficients 
which are bounded functions, so it is a lipschitzian function on bounded subsets 
and satisfies lemma hypothesis. □ 

2.2 Forced vibration, double scale expansion of order 1 
2.2.1 Derivation of the expansion 

Here we consider a similar system with a sinusoidal forcing at a frequency close 
to the free frequency (so called primary resonance); in the linear case, without 
damping, it is well known that the solution is no longer bounded when the 
forcing frequency goes to the free frequency. Here, we consider the mechanical 
system of previous section but with periodic forcing and we include some light 
damping term; the scaling of the forcing term is chosen so that the expansion 
works properly; this is a known difficulty, for example see [31] . 

du"^ 

ii + uTu + eXii + cu 2 H — = e 2 F cos(uj f t). (24) 

We assume positive damping, A > and excitation frequency (I> e is close to an 
eigenfrequency of the linear system in the following way: 

d> e = oj + ea. (25) 

We look for a small solution with a double scale expansion; to simplify the 
computations, the fast scale T is chosen e dependent and we set: 

T = u> e t, Ti = et and D u = Dm = J^-, (26) 

01 o Oil 

so 

— = uiD u + eDm and — - = oj 2 D 2 u + 2eu} e D Q D 1 u + e 2 D\u- 1 (27) 
at at 2 

equation (|2"5"j) provides 

a; 2 = lu 2 + 2eujo- + e 2 a 2 . (28) 
With ((23 (25), (@7|, ® and the ansatz 

u = eu 1 (T ,T 1 ) + e 2 r(To,T 1 ,e), (29) 
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we obtain: 

du du\ 9 dr dui , „ 9/ dr . , . 

Tt=^ + f dt= f HT + *" Dor + e ( dT - = ( 3 °) 

elCjDom + eDiUi] + e 2 ujD r + e 2 (— - cjZV) = (31) 

efwDoUi + €(tD ui + eD\U\\ + ( 2 u)D$r + e 2 (- u;D r), (32) 

di 

where we remark that ^ — u)D$r = eaD^r + eD\r is of degree 1 with respect 
to e. For the second derivative, as for the case without forcing, we introduce 

1 fd 2 r \ 
V 2 r = - ( —j — urD^rj , with the expansion (33) 

V 2 r = 2oj[aD 2 r + D D ir ] + e [a 2 Dlr + 2aD D ir + D\r] , (34) 



d 2 u d 2 ui ,d 2 r d 2 i*i 99^9 q„ 

d^ = ^ + e V - ^ + *" D ° r + eV2r (35) 

= e [£j 2 D 2 Ui + 2euD Q D lUl + € 2 D\ Ul ] (36) 

+ e 2 u 2 D 2 Q r + e 3 V 2 r (37) 

= e{uj 2 D 2 Ul + 2euo (aD^Ut + D D lUl ) + (38) 

e 2 [a 2 Dim + 2aD D lUl + D\ Ul ] } (39) 

+ e 2 u 2 Dlr + e 3 V 2 r; (40) 

the last term in the right hand side will be part of the remainder R of 
equation ((42j . We plug previous expansions into ([24)) ; we obtain: 

f u 2 (D 2 ,u 1 +ui) = 0, 
1 £)2 r + r=§, with 



5 2 = - {cu 2 + du\ + 2u[D Q D 1 ui + cr Dim] + XuiDqUx} (42) 
+fcos(To)-eiJ(«i,r,c) (43) 

and with 

R(ui, r, e) = D\u\ + 2cu%r + idu\r + <r 2 

£>gui + 2aI>o£>iUi+ (44) 
\(wD r + (tDqUx + D1U1) + D 2 r (45) 
dr 

+X{—-ujD Q r) + ep(u 1 ,r,e). (46) 

Set 0(T o ,Ti) = T + We solve the first equation of gTJ : 

ui = a(Ti)cos((9) (47) 
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then we use T = 6(T , T x ) - /3(7i), and we obtain 



S 2 = -^-(1 + cos(20)) - — (cos(30) + 3cos(0)) + 

2uj(Dia s'm(9) + aD\f5 cos(9)) + Ioujo, cos(#) + a\uj s\n(9) 

+ F sin(0) sin(/3(Ti)) + F cos(0) cos(/3(Ti)) - eR(«i, r, e) 



(48) 



or 



5 2 = [2wDia + Aaw + FsmfJ))] sin(0) 



2ujaDif3 + 2craja 



3da J 



Fcos(/3) 



cos(6>) 

- eR(m,r, e) 



with 



S| = -^(1 + cos(20)) - ^ (cos(30)) ; 



(49) 
(50) 



note that S\ is a periodic function with frequency strictly multiple of 1. 



Orientation. By enforcing 



2u)D\a + Xauj = —Fsm(f3), 
2ciluDi/3 + 2<jLua 



Ma' 



-Fcos(/3), 



(51) 



S*2 = S\ — eR(u\, r, e) contains neither term at frequency 1 nor at a frequency 
which goes to 1; this point will enable to justify this expansion under some 
conditions; first, we study stationary solution of this system and the stability of 
the dynamic solution in a neighborhood of the stationary solution. 



2.2.2 Stationary solution and stability 

Let us consider the stationary solution of ([5Tj) . it satisfies: 



UXau + FmaiP)] =0, 



2, 
1 

2u> 



2>jjo 



3da^ 
4 



Fcos(ff) 



= o, 



(52) 



Now, we study the stability of the solution of ([51~T) . in a neighborhood of this 
stationary solution noted (a, j3); set a = a + a, j3 — ft + j8, the linearized system 
is written 



J 



manipulating, we obtain the jacobian matrix 
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The matrix trace is —A, and the determinant is 

, /T . A 2 ,9da 2 w 3da 2 
det(J) = — - (— - a)(a 



4 v 8a; /v 8w " 

we notice that the determinant is strictly positive for a = so by continuity, it 
remains positive for a small; moreover ■4-det(J) < for <r < so det(J) > 
for er < 0; by studying the trinomial in a, we notice that the determinant is 

positive when this semi-implicit inequality is satisfied: a < — \\J 9 iqS ~ 
so in these conditions, the two eigenvalues are negative; then the solution of the 
linearized system goes to zero; with the theorem of Poincare-Lyapunov (look 
in the appendix for the theorem 15.11 ) when the initial data is close enough to 
the stationary solution, the solution of the system (|51|) . goes to the stationary 
solution. We expand this point, set 

-Xau -Fsin(/3) \ 
G{y)=\ _ f ,.,„ _ M*L\ Fcosjp) ); (54) 



PJ \ - 2wcr- 



the system (I5^|) may be written y = G(y); denote y = (|), the solution of (|52j): 
perform the change of variable y = y + x, we have G{y +x) = G(y) + Jx + g(x), 
with g{x) — o(||x||); the theorem 15. II may be applied with A — J, B — 0, here 
the function g does not depends on time. 



Proposition 2.2. If a < -1^/2^ _ \2 , the stationary solution of lbl\) 
is stable in the sense of Lyapunov (if the dynamic solution starts close to the 
stationary solution of $52\) . it remains close to it and converges to it ); to the 
stationary case corresponds the approximate solution of \2J$ U\ — a cos(To + j3) , 
it is periodic; for an initial data close enough to this stationary solution, u\ = 
a(Ti)cos(To + /3(Ti)), with a, j3 solutions of i51\) : it goes to the solution h52\) 
a,/3 when T\ — > +oo. 

With this result of stability, we can state precisely the approximation of the 
solution of (|24|) by the function u\. 



2.2.3 Convergence of the expansion 
Proposition 2.3. Consider the solution of \21$ with 

u(0) = ea + o(e), u(0) = -ecoa sin(/3 ) + o(e), 
with ao,/3o close of the stationary solution (a,/3), 

|ao-o| <eC u \fa-P\ <eC 2 . 

3da 2 

the following expansion is satisfied 

u(t) = ea(et) cos(w e i + (3(et)) + e 2 r(e,t), 



}}-!,< ii n < 2gl _ I^sgi _ A2, there exists 7 > such that for all t<t e = Z, 
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with ui e — uj + ecr and r uniformly bounded in C 2 (0,t e ) and with a, /3 solution 

of m>- 



Proof. Indeed after eliminating terms at frequency 1, we go back to the variable 
t for the second equation (jITjl 

S 

r + uj 2 r~^r, with (55) 
S 2 = Si[t,e)-eR{u x ,r,e), (56) 

with 

dr 

R(u u r, e) = R(ui,r, e) - V 2 r - A(— - uD r), (57) 
with all the terms expressed with the variable i; we have 

S«(t, e ) = ~ ca2 ( et ) (i + cos(2(a; e t + yg(et))) - ^L^£i ( CO s(3(a) 6 t + /?(e*)) , (58) 

this function is not periodic but is dose to the periodic function: 

S\{t, e) = ^-(1 + cos(2(w e i + p)) - ^- (cos(3(w e i + 0)) (59) 



and for i < 2 as the solution of ([5~Tj) is stable: it remains close to the 
stationary solution 

\a(et) - a\ < eCx, |/3(et) - p\ < eC 2 (60) 

and 

\Si-S\\<eCz- (61) 

so this difference may be included in the remainder R. We use lemma 15.11 with 
S = S\\ it satisfies lemma hypothesis; similarly, we use g = R; it satisfies the 
hypothesis because it is a polynomial in the variables r,U\,e, with coefficients 
which are bounded functions, so it is lipschitzian on bounded subsets. □ 

Remark 2.1. Under different hypothesis and for systems involving first order 
derivatives, a result of convergence may be found in !32\j . 

2.2.4 Maximum of the stationary solution, primary resonance 

Consider the stationary solution of (I5T1) . it satisfies 



Xatu = -Fsin(/3), 
a[2ojcr-^f) =-Fcos(f3), 



(62) 



manipulating, we get that a is solution of the equation: 



/(«.,, ] =\2 a 2 u 2 + a 2 (2u J o--^-) 2 -F 2 =0. (63) 
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We compute 



& f 

— — = Aa 2 uj(2uxj 
da 



Ma 2 



d f o \2 2 
— = 2aX uj 

da 



2ii i 2ujo — 
da 2 



Zda 2 



da 3 



3da 2 



ia 2 uj 2 . 



(64) 

— ()—— ( 2^'rr — — — | . (65) 

(66) 
(67) 

For a close enough to the solution of l|j = 0, ^ is small, S is not zero, and 
with the implicit function theorem this equation defines a function a(a); lets 
use : 

8a U- , 9 2 a 
5a " # and a(j2 



a<r 2 
a/ 

(9 a 



In our case, when 4p = 0, we have 



3da 2 
8w 



9/ 
9<z 



= 2aA 2 w 2 , 



<9ct 2 



= 8a 2 cA 



(68) 



so the second derivative |g§ < and a is maximum at the frequency of the free 
periodic solution. 

Proposition 2.4. The stationary solution of \51\) satisfies 



2aL>jo 



Xacu +Fsin(/3) : 
^ +Fcos(/3) 



0, 

= 0, 



it reaches its maximum amplitude for a = and /3 = § 
is at the angular frequency 



(69) 

kit; the excitation 



3e- 



d^_ 

&7 



id F = Xtua: 



it is the angular frequency v e of the free periodic solution h23\) for this frequency, 
the approximation ( of the solution up to the order e) is periodic: 



p 

u(t) = e— sin(cD e i) + e 2 r(e, t). 



(70) 



Xui 

Remark 2.2. We remark that this value of a = 3 4 2l - is indeed smaller than the 
maximal value that a may reach in order that the previous expansion converges 
as indicated in proposition [ 



Remark 2.3. We note also that when the stationary solution reaches its max- 
imum amplitude we have F — Xuia and so we can recover the damping ratio X 
from such a forced vibration experiment; this is a close link with the linear case 
(see for example JSSjj or the English translation JS^j). This is quite interesting 
in practice as the damping ratio is usually difficult to measure; we have here a 
kind of stability result for this experiment. 
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3 System with strong local cubic non linearity 

In the previous section, we have derived a double scale expansion of a solution 
of a one degree of freedom free vibrations system and damped vibrations with 
sinusoidal forcing with frequency close to free vibration frequency. Now, we 
extend the results to the case of multiple degrees of freedom. 

3.1 Free vibrations, double scale expansion 

We consider a system of vibrating masses attached to springs: 



The mass matrix M and the rigidity matrix K are assumed to be symmetric 
and positive definite. We assume that the non linearity is local, all components 
are zero except for two components p — 1, p which correspond to the endpoints 
of some spring assumed to be non linear: 

$ p -i(u, e) = c(u p - Up_i) 2 + ^(u p - Up_i) 3 , 5>p = -$ p _i, p = 2,...,n. (72) 

If the non linear spring would have been the first or the last one, the expression 
of the function $ would depend on the boundary condition; each case would 
be solved using the same method with slight changes in some formulas. In 
order to get an approximate solution, we are going to write it in the generalized 
eigenvector basis: 



Mu + Ku + $(u, e) =0. 



(71) 



Ktfik = LUkM(j)k, with 4>\M(f>i — Ski, k, I = 1 . . . , n. 
So we perform the change of function 



(73) 



n 




(74) 



we obtain 



n 




k = 1 . . . , n. 



(75) 



i=i 



As $ has only 2 components which are not zero, it can be written 



n 




k = 1 . . . ,n, 



(76) 



or more precisely 





= 0, k = l...,n. (77) 
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As for the 1 d.o.f. case, we use a double scale expansion to compute an approx- 
imate small solution; more precisely, we look for a solution close to the normal 
mode of the associated linear system; we denote this mode by subscript 1; ob- 
viously by permuting the coordinates, this subscript could be anyone (different 
of p, this case would give similar results with slightly different formulas); we set 



and we use the ansatz. 



so that 

d 2 y k 
dt 2 

with 



T = u>it, T x = et (78) 
y k = ey\ (T , Tj ) + e 2 r k (T , Tj , e) , (79) 

eu 2 D 2 yl + e 2 [2Lu 1 D D 1 y 1 k + wlD 2 r k ] + e 3 [T^yi + V 2 r k ], (80) 



V 2 r k = 1 C^- - Lj 2 D 2 r^J = 2ui 1 D D 1 r k + (D 2 r k . (81) 

We plug previous expansions into (|77[) . By identifying the coefficients of the 
powers of e in the expansion of (1761) . we get: 



J <4Dly\+u\yl = 0, k= 1 . . . , n, 

1 uj\D 2 ru + ^\xk = S 2 ,k, k = 1 . . . ,n, with 

to simplify, the manipulations, we set 

8p<f>l = {4>l,p ~ 0Z,p-l), 

so: 

-8 p (j) k 



(82) 



S2, k 

with 
and 

S2,k 



*p-i + ^ r ^^j ~ ^iD D iy l - eK k , (83) 



TZ k = {D 2 yl + V 2 r k ) 



(84) 



-5 p 4> k 



-2uj l D Q D l y 1 k -en k . (85) 



The formula may be expanded 



S 2 ,k = -8: 



• P (p k 



- 2^ x D^D iy \ -eR k (y 1 ,r,e), (86) 
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where 
R k (y\r,e)=K k 



+ S p (/>k 



i,3 



{^ViVj r l + 3eylr 3 ri + S^nrjri) SpfaSpfijSpfii 

ijl 



(87) 



We set 0(T o ,Ti) = T + /3(T X ) and we note that D o = 1, L>i0 = £>i/9; we 
solve the first set of equations ([52)1 , imposing O(e) initial Cauchy data for k ^ 1; 
we get: 

= a (Ti) cos(0), and y^ = 0(e), k = 2 . . .n; (88) 
we put terms involving y^, fc > 2 into i?^; so we obtain 



S2,i=-5 p (j)i c (ylSpcj)!) 2 + d (ylSp(t)iy 



2 Wl A) Ayi - eiZi (y 1 , r, e) and (89) 



S" 2: fc = -£ p <fo c (yi<5 P 0i) 2 + d (ylSpcf)!) 3 



Using (|55]l. we get: 



ei? fe (y 1 ,r,e)forfc^l. (90) 



Sa,i = -<5 P ^i [^i(l + cos(20)) (5^x) 2 + 



da\ 



((cos(30) + 3cos(0))(<5 p </>i) 3 ) 



2wi(Diai sin(0) + a x £)i/3i cos(0)) - eiJ^y 1 , r, e) and (91) 



S 2 ,fc = -<5 P 0fc [^(1 + cos(20)) (V^) 2 + 



da? 



((cos(30) + 3cos(0))(5 p 0i) 3 



-ei? fc (y 1 ,r,e)forfc^l. (92) 
We gather the terms at angular frequency 1 in 52,1 : 

■da\. 



#2,1 = S> 



-3cos(0)(5 p (^i)' : 

+ 2wi(Diai sin(0) + ai-Dift cos(0)) + - ei?(y\ r, e), (93) 
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with 

°2,1 — 

If we enforce 



-S p 4>i 



^-(1 + cos(20))(^0 1 ) 2 + cos^X^) 3 



(94) 



.Didi = 0, and 2uJxaxDiPx = {5 P 4>\) 



4 3da 3 

— , so that 



the right hand side 



ai = ax.o, fix = Px,oTx with /3i 



5 2 ,i =S{ 1 -eR 1 (y 1 ,r,e) 



3da 2 

~&7 



(Mi) 4 Ti, (95) 



(96) 



contains no term at angular frequency 1; for the other components, without any 
manipulation, there is no trouble with the frequencies if we assume that all the 
eigenfrequencies uj k for k = 2 . . . n are not multiple of ujx (oJk ^ Q^i for q = 1 or 
Q '/ 

In order to prove that r is bounded, after the elimination of terms at fre- 
quency 1, we write back the equations with the variable t, for the second set of 
equations of I 



with 
where 

S 2,l = -5p<t>i 

and 

S2M = -5 P 4>k 



ufr k + wlr k = §2,k, for k = 1, ...n, 
§2,1 = Si x - e-Ri(j/V.e): 



-^(1 + cos(2(^ + p lfi et)){d p ^) 2 



(97) 
(98) 



da 3 

+ -ji cos(3( Wl < + Px.^t)){8 p 4>xf 



(99) 



l Ml + cos(2( Wl t + 0i, o rf))) (Mi)' + 



da 3 

-^i ((cos(3(wit + ft l0 ct)) + 3 cos((wi* + /3 ljO et)))(<5 P 0i) 3 ) 



and where 



-e J R fe ( 2 / 1 ,r,e)forfc^l (100) 
R k (y\r, e) = R k {y\r, e) - V 2 r k . (101) 
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Proposition 3.1. Under the assumption that Uk andu\ are Z independent for 
k =/= 1, there exists 7 > such that for all t < t e — % } the solution of |76| ) with 
initial data 

lft(0) = coi, 0) tfi(0) = 0, y fc (0) = O(e 2 ), fo(0) = (102) 
satisfy the following expansion 

da 2 

yi(t) = eaocos(v e t) + e 2 n(e,t) with v e = oj\ + 3e— —(<pi, p ~ 0i. p -i) 4 , (103) 

2/ fe (t) = e 2 r fc (e,t), (104) 

wii/i r,t uniformly bounded in C 2 (0, i e ) /or k = 1, . . . n and u)i, <pi are the eigen- 
value and eigenvectors defined in J7ff| ). 

Corollary 3.1. The solution of {77} (7J§) with 

^«(0) = eai,o, <^(0)=0, <^u(0) = 0(e 2 ), ^u(O) = 0, 
zyzi/i LJk,<Pk are the eigenvalue and eigenvectors defined in (73j) 

n 

is «(t) =J2vk(t)<pk, (105) 
fc=i 

£/ie expansion of of previous proposition. 

Proof. For the proposition, we use lemma [5T4l Set Si = S| Sfc = S2,fc for 
fc = 1, . . . n; as we have enforced (|95|) . the functions <Sfc are periodic, bounded, 
and are orthogonal to e ±lt , we have assumed that uik and wi are Z independent 
for fc 7^ 1; so Sk, k = l,...,n, satisfies the lemma hypothesis. Similarly, set 
o = i?, its components are polynomials in r with coefficients which are bounded 
functions, so it is lipschitzian on the bounded subsets of R, it satisfies the 
hypothesis of lemma 15.41 and so the proposition is proved. The corollary is an 
easy consequence of the proposition and the change of function (|107[) □ 

Remark 3.1. We have obtained a periodic asymptotic expansion of a solution 
of system J7i[ ), |7£[ ); they are called non linear normal modes in the mechanical 
community ([23, 221. In the next section, we shall derive that the frequencies of 
the normal mode are resonant frequencies for an associated forced system, the 
so called primary resonance; secondary resonance could be derived along similar 
lines. 

3.2 Forced, damped vibrations, double scale expansion 
3.2.1 Derivation of the expansion 

We consider a similar system of forced vibrating masses attached to springs with 
a light damping: 

Mil + eCu + Ku + e) = e 2 Fcos(w £ i), (106) 
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with the same assumptions as in subsection 13. II We assume that the frequency 
of the driving force is close to some frequency of the linearised system (primary 
resonance); we denote this frequency with the subscript 1: oj e — + ea 

We assume that the non linearity is local, all components are zero except 
for two components p — 1, p which correspond to the endpoints of some spring 
assumed to be non linear. As for free vibrations, we perform the change of 
function 



k=l 



(107) 



with 4>k , the generalised eigenvectors of (|73|) . As the damping matrix C is usually 
not well defined, to simplify, we assume that it is diagonal in the eigenvector 
basis 4>k, k = 1, . . . n. We obtain 



Vk + (^kVk +uiyk + 4>l^(^2yi(j)i,e) = e 2 f k cos(Q e t), k = l...,n, (108) 



1=1 



with fk — 4>kF. As for the free vibration case, $ has only 2 components which 
are not zero, so the system can be written 

n 

Vk + eA fe y fc +0JkVk + {<t>k,p-i - <l>k,p,e) ^p-ij/^Vi^i) = e 2 f k cos(u t t), 



k = l...,n, (109) 



or more precisely 



yk + e\kVk + ulllk + {4>k, P -i - 4>k, P ) 



y^ofc.p - <t>j,p-i) j + 

\»=1 / 



yi(<Pi,p - Pi,p-i 



e 2 f k cos(io e t), 



k = l...,n. (110) 



As for the 1 d.o.f. case, we use a double scale expansion to compute an approx- 
imate small solution; we use a fast scale which is e dependent; we set 



T = & e t, Ti = et, 

and we use the "ansatz" 

Vk = «/fc(T ,Ti) + e 2 r k (T ,T 1 ,e), 

so that 



(111) 
(112) 



dy, 



dr k 



-j^ = e^Doyl + eaDoyl + eD l y 1 k ] + e 2 u x D r k + f 2 {-£ - ^D r k ), (113) 



dt 



3 SYSTEM WITH STRONG LOCAL CUBIC NON LINEARITY 18 



(l/J = eUiDlyl + 2eu! [aD 2 yl + D Diyl] + 

e 2 [a 2 D 2 yi + 2*D D lV l + D 2 yl] J 

+ e 2 cu 2 1 D 2 r k +e 3 V 2 r k , (114) 



with 



1 f d 2 r k 
e \ di 2 



V ^ = " ( - ^iD 2 r k ) = 2u>i(aD 2 r k + D D x r k ) 



e [a 2 D 2 r k + 2aD D irk + D 2 r k ] . (115) 



We plug previous expansions into (|110l) . By identifying the coefficients of the 
powers of e in the expansion of (|110l) . we get: 



u} 2 D 2 yl+u 2 k yl=0, k=l...,n, 

u\Dlr k + io\r k = S 2 , k , k = 1 . . . , n, with 



(116) 



5; 



a,* = (^(eyl + e 2 n)^,ej+2Lj 1 [D D 1 yl+aD 2 yl}+X k Lj 1 D yl 

+ .f k cos(T )-eR k (y 1 ,r,e), (117) 



where we gather higher order terms in R k and to simplify, the manipulations, 
we have set 

5 p 4>i = (<f>i tP - 4>i, p -i), 



so: 



'2, k 



Sp4> k 



- 2u 1 [D D 1 y l k + aD 2 yl] - X^DqvI 

+ f k cos(T )-eR k (y 1 ,r,e). (118) 



The formula may be expanded 



S 2 ,k = -5 P 4>k 



'^r,yly)Sp<i>iSp<i>j + d^y}y)y} SpfaSpfaSpC, 

- 2cj 1 [D D l yl + aD 2 yl] - \ k ^D Q yl 

+ f k cos(T )-eR k (y\r,e). (119) 
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We set 0(T o ,Ti) = T + £(T X ) and we note that A)0 = 1, D X B = Dx/3; we 
solve the first set of equations (jl 16[) . imposing initial Cauchy data for k ^ 1 of 
order O(e) we get: 



vl 



ai(Ti)cos(6»), and = O(e), k = 2...n 
we put terms involving y\ into Rk for A: > 2 and so we obtain 



(120) 



£2,1 = S. 



2wit.Do.Diyi 1 +<?D 2 yl} - A lW iA>l/i + A cos(T ) - eR^^e) and (121) 



S2,k = —S p (j)k 



c(ylS P (f>i) 2 +d{y\5 p (j)i 



f k cos(To) - eR k (y\r, e) for k ± 1. (122) 



Using (TT2171) . we S et: 



52,1 = -<5 P 0i [-^(1 + cos(20)) {Sp^f + 



((cos(36>) + 3cos(0))(<5 p< />i) 3 ) 



2wi[D>iai sin(0) + aiDifii cos(8) + aa x cos(#)] + AiWidi sin(#) 

+ /i[sin(0)sin(/3) + cos(6») cos(/3)] - ei^y 1 , r, e) and (123) 



5 2 , fe = ~5 P 4>k [^(1 + cos(20)) (<5 p 0i) 2 + 
da 3 

-ji ((cos(36») +3cos(6»))((5 p 0i) 3 ) 
+ / fe [sin(0) sin(/3) + cos(6>) cos(/3)] - eR k (y\r, e) for fc ^ 1. (124) 

We gather the terms at angular frequency 1 in 52.1 



52,1 = 5p4>l 

with 
at 

°2,1 — 



da 3 

~ cos(6»)(^0i) 3 + 2io 1 (a 1 D 1 (3 1 + aai) + fx cos(/3) 



cos(6>) 



wi(2Dio x + Aiai) + A sin(/3)l sin(0) + S\ 1 - eR(y\r, e), (125) 



-5 p <j>x 



CO 



da\ 



^(1 + cos(20))(a„0i) 2 + cos(30)(5 p 0i; 



(126) 
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Orientation If we enforce 



ui 1 (2D 1 ai + Aiai) = -fx sin(/3i), and 
2u 1 (a 1 D 1 l3 1 + aax) = ^{8 p (j)i) 4 - /icos^i), 1 j 

the right hand side 

£u = 4i- ei? i(y 1 ' r ' e ) ( 128 ) 

contains no term at angular frequency 1; for the other components, without any 
manipulation, there is not such terms , if we assume that all the eigenfrequencies 
ujk for k — 2...n are not multiple of lj\ (ui k ^ qtui for q = 1 or q = 2, 
q = 3). This will enable us to justify this expansion; previously, we study the 
stationary solution of this approximate system and the stability of the solution 
in a neighbourhood of this stationary solution. 



3.2.2 Stationary solution and stability 

The situation is very close to the 1 d.o.f. case; except the replacement of d by 
of d = d(5 p (j)i) 4 , the system (|127p is the same as (|5T|) ; the other components are 
zero. We state a similar proposition 

Proposition 3.2. When a < — \\J ^iqS ~ ^i > ^ e s ^ a ^ onar y solution of 
1127\ ) is stable in the sense of Lyapunov (if the dynamic solution starts close 
to the stationary one, it remains close and converges to it); to the stationary 
case corresponds the approximate solution of |77| y\ — a i cos (To + y\ = 

0(e), k = 2, . . . ,71, it is periodic; for an initial data close enough to the sta- 
tionary solution, y\ = a{T\) cos(To + @i{Ti)), y\ — 0(e), k — 2, . . . , n with 
a, /?i solutions of \12 r /\ ) with d replaced by d ; they converge to the stationary 
solution ai,f3i when T\ — > +oo. 



3.2.3 Convergence of the expansion 



In order to prove that r is bounded, after the elimination of terms at frequency 
1, we write back the equations with the variable t, for the second set of equations 
of (JH2D- 

■>.. ■> ,~, 

.n, 



with 
where 

S 2.1 = -tpfa 



u\r k + u? k r k = S 2 ,k for fc = 1, 
5*2,1 = 5*2,1 - e -Ri(yV, £ ): 



c(ai(et)) 2 



(1 + cos(2(u! e t + 0i(et))(6 p (f>i 
da\ 



(129) 
(130) 



+ _i CO s(3(w £ i + /3i(ei))(^0i) c 



(131) 
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and 

S2,k = —5 p 4>k 



c{ai< f t)2 (1 + cos(2(c^ + p 1 (et))) (S^j 2 - 



da 3 

-± ((cos(3(w e i + ^ ( et)) +3cos((w £ t + /3 1 (et)))(<5 p 1 ) 3 

-eR k {y 1 ,r,e)iovk^l, (132) 



where 



^fc(yV, e) = i? /c (y 1 ,r, e) - P 2 rfc - -M^r ^ u k D r k )- (133) 



Proposition 3.3. Under the assumption that uj k and are Z independent for 
k ^ 1, t/iere exists 7 > suc/i £/ia£ /or all t < t e = ^ , the solution of illO\) with 
initial data 

yi(0) =eai, + O(e 2 ), ^(0) = -ewai.o sin(/3 1)0 ) + 0(e 2 ), (134) 
y k (0) = O(e 2 ), y k (0)=0 (135) 
and with the initial data close to the stationary solution 
|ai,o - ai| < eCi, |/3i )0 - < eCi 
satisfy the following expansion 

yi(t) = eai(et) cos(w e t + /3i(et)) + e 2 ri(e, i), wifft (136) 
Mk(t) = e 2 r-fe(e,i), (137) 

wii/i ai,/3i solution of \121^ and with r k uniformly bounded in C 2 (0,t e ) for 
k = 1, ...n and u)y,<j>\ are the eigenvalue and eigenvectors defined in \73\ and 
a%Pi are solution of \12 r ty 

Corollary 3.2. The solution of (TUS\) . (7|) with 

i#u(0) = eai.o, 0f«(O) = -ewio 1 , O 8inG8i,o), $u(0) = 0(e 2 ), <^u(0) = 0, 
with ujk^k the eigenvalues and eigenvectors defined in [73jj , 

n 

is u(t) =J2vk{t)^k, (138) 
fe=i 

with the expansion of y k of previous proposition. 

Proof. For the proposition, we use lemma I5~4l Set Si = S| u Sfc = S2, k f° r & = 
1, . . . n; as we have enforced (|95p. the functions Sjfe are periodic, bounded, and are 
orthogonal to e ±4 *, we have assumed that uj k and uj\ are Z independent for k ^ 1; 
so 5 satisfies the lemma hypothesis. Similarly set 5 = _R, it is a polynomial 
in r with coefficients which are bounded functions , so it is lipschitzian on 
the bounded subsets of R, it satisfies the hypothesis of lemma 15.41 and so the 
proposition is proved. The corollary is an easy consequence of the proposition 
and the change of function (|107|) □ 
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3.2.4 Maximum of the stationary solution 

As equation (|127|) is similar to the equation (|5Tj) of the 1 d.o.f. case, we get also 
that the stationary solution reaches its maximum amplitude to the frequency of 
the free periodic solution. 

Consider the stationary solution of (|127p . it satisfies 

f Aiai^i = -/isin(/3i), 

{ afaa-^f) =-/ 1 cos(/3 1 ); (139) 
manipulating, we get that a\ is solution of the equation: 

f{a 1 ,a) = \la\^ + alLu l a-*-^ - f 2 = 0. (140) 
As for the 1 d.o.f case, we can state: 

Proposition 3.4. The stationary solution of {127^ satisfies 

J Aiaicji +/i sin(/3i) = 0, 

) 2aiWl a-^ +/ 1 cos(/3 1 ) = 0, 



(141) 



it reaches its maximum amplitude for a = ^jf- an d ^ — TL _j_ J^^- tfo e excitation 
is at the frequency 

da 2 

Q e = u>i + 3e — -, with d = d(<f>i tP — 3>i )P _i) and F — AiWiai 

where $i is the eigenvector of the underlying linear system associated to uj\; 
uj e is the frequency of the free periodic solution \23\) : for this frequency, the 
approximation ( of the solution up to the order e) is periodic: 

r 

yi(t) = e-^sm(u,t)+e 2 r(e,t), (142) 
Aicji 

y k (t) = e 2 r k (e,t). (143) 
As for the 1 d.o.f. case we can remark the following points. 

Remark 3.2. This value of a = ^h. j s indeed smaller than the maximal value 
that a may reach in order that the system be stable and that the previous expan- 
sion converges as indicated in proposition \2.3l 

Remark 3.3. We note also that when the stationary solution reaches its max- 
imum amplitude we have f± = AiWiai and so we can recover the damping ratio 
Xi from such a forced vibration experiment; this is a close link with the lin- 
ear case (see for example or the English translation 134-1)- This is quite 
interesting in practice as the damping ratio is usually difficult to measure. Ob- 
viously, we can recover the damping ratio for other frequencies by performing 
other experiments. 

We can also consider this result as a stability of the process used in the linear 
case with respect to the appearance of a small non-linearity. 
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4 Conclusion 

For differential systems modeling spring-masses vibrations with non linear springs, 
we have derived and rigorously proved an asymptotic approximation of periodic 
solution of free vibrations (so called non linear normal modes); for damped vi- 
brations with periodic forcing with frequency close to free vibration frequency 
( the so called primary resonance case), we have obtained an asymptotic ex- 
pansion and derived that the amplitude is maximal at the frequency of the non 
linear normal mode. 

Acknowledgment We thank S. Junca for his stimulating interest. 

5 Appendix 

Lemma 5.1. Let w e be solution of 

w" +w = S(t,e) + eg(t,w,e), (144) 
io(0)=0, to'(0)=0. (145) 

If the right hand side satisfies the following conditions 

1. S is a sum of periodic bounded functions: 

(a) for all t and for all e small enough, S(t, e) < M, 

(b) Jq™ e lt S(t,e)dt = 0, e~ lt S(t,e)dt = uniformly for e small 
enough 

2. for all R > 0, there exists such that for \u\ < R and \v\ < R, the 
inequality \g(t,u,e) — g(t,v,e)\ < kji\u — v\ holds and \g(t,0, e)\is bounded; 
in other words g is locally lipschitzian with respect to u. 

then, there exists 7 > such that for e small enough, w e is uniformly bounded 
in C 2 (0,T £ ) with T e = 2 

Proof. The proof is close to the proof of lemma 6.3 of [TS]; but it is technically 
simpler since here we assume g to be locally lipschitzian with respect to u 
whereas it is only bounded in |15) . 

1. We first consider 

Wi" +wi =S(t,e), (146) 
Wl (0)=0, w[(0) = 0: (147) 

as S is a sum of periodic functions which are uniformly orthogonal to e lt 
and e~ lt , wi is bounded in C 2 (0,+oo). 
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2. Then we perform a change of function: w = w\ + w 2 , the following equal- 
ities hold 



w 2 v + w 2 = eg 2 (t,w 2 ,e), 
w 2 (0) = 0, w' 2 (0)=Q, 



(148) 
(149) 



with g 2 which satisfies the same hypothesis as g: 

for all R > 0, there exists kpi such that for \u\ < R and \v\ < R, the follow- 
ing inequality holds \g 2 (t,u,e) — g 2 (t, v,e)\ < kn\u — v\. Using Duhamel 
principle, the solution of this equation satisfies: 



w 2 = el sin(f — s)g 2 (s,w 2 (s), e)ds 



(150) 



from which 



\w 2 (t)\<e \g 2 (s,w 2 (s),e)-g 2 (s,0,e)\ds + e \g 2 (s, 0, e)|ds (151) 
Jo Jo 

so if | iu | < R, hypothesis of lemma imply 

\w 2 (t)\<e[ k R \w 2 \ds + eCt. (152) 
Jo 

A corollary of lemma of Bellman-Gronwall, see below, will enable to con- 
clude. It yields 

M*)| < r- (exp(eM) - 1) ■ (153) 
Now set T c — swp{t\\w\ < R}, then we have 

C 

R<T~ (exp(ek R t) - 1) ; 

this shows that there exists 7 such that \w 2 \ < R for t <T £ , which means 
that it is in L°°(0,T e ) for T e = ^; also, we have w in C(0,T e ) then as w is 
solution of (|144j) . it is also bounded in C 2 (0,T £ ). 

□ 

Lemma 5.2. (Bellman-Gronwall, [35, 361) Letu,e,/3 be continuous functions 
with P>0, 



u(t) < e(t) + / /3(s)u(s)ds forQ<t<T, 



then 



u(t) < e(i) + / /3(s)e(s) 



exp( / P{t)cIt 



ds. 



(154) 
(155) 
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Lemma 5.3. ( a consequence of previous lemma, suited for expansions, see 
[121) Let u be a positive function, 82 > 7 8\ > and 

u(t) < S 2 t + Si ( u(s)ds, 
Jo 

then 

u{t) < ^ (exp(ftt) - 1.) 
Lemma 5.4. Let v e = [v{, . . . , v^] T be the solution of the following system: 

<4W +"tv e k = S k (t) + eg k (t,v e ). (156) 

If oji and u>k o,re Z independent for all k — 2 . . . N and the right hand side 
satisfies the following conditions with M > 0, C > prescribed constants: 

1. Sk is a sum of bounded periodic functions, \Sk{t)\ < M which satisfy the 
non resonance conditions: 

2. Si is orthogonal to e ±lt , i.e. f Q ™ Si(t)e ±lt dt — uniformly for e going to 
zero; 

3. for all R > there exists kn such that for \\u\\ < R, \\v\\ < R, the following 
inequality holds for k = 1, . . . , N : 

\g k (t,u,e) - g k (t, v,e)\ < k R \\u - v\\ 

and \gk(t,0,e)\ is bounded 

then there exists 7 > such that for e small enough v e is bounded in C 2 (0,T e ) 
with T e = Z 

Proof. 1. We first consider the linear system 

wf(ufe,i)" +u>lv kA = S k , (157) 
u M (0)=0and(« M )' = 0. (158) 

For k = 1, with hypothesis l.a, Si is a sum of bounded periodic functions; 
it is orthogonal to e ±%t , there is no resonance. For k 7^ 1, there is no 
resonance as ^ ^ Z with hypothesis l.b. 

So v k ,i belongs to C^ & '(/,%) for k — l,.,.,n. 
2. Then we perform a change of function 

V% = V k ,l + Vfc, 2 

and v e k 2 are solutions of the following system : 
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^i(vk,sT + ^t v k,2 = egk,2{t, vk,2, e), k = l,...,N, (159) 
«fc,a(0) = 0, («fc l2 )' = 0, fc = l,...,iV, (160) 

with 

where g^.2 satisfies the same hypothesis as gk- 

for all R > there exists Ur such that for || Uk ||< -R, || ffc ||< -R, the 
following inequality holds for k = 1, . . . , N : 

II 9k,2(t,u k ,e) - gk,2(t,Vk,e) \\< k R \\ u k - v k \\ . (161) 
Using Duhamel principle, the solution or the equation (|159[) satisfies: 

v t,2 = e [ sin (*- s )5fc,2(s,< 2 (s),e)ds, (162) 
Jo 

so 

II v*,2(*) ll< « / lls*,2(a,«fc l2 (*),e)-5fc,2(a.0.e)||da+ 
Jo 

e f \\g kt2 (s,0,e) \\ ds, (163) 
Jo 

so with (|161[) , we obtain 

Ka(*) ll<£ f k\\vl 2 {t) \\& s + eCt. (164) 
Jo 

We shall conclude using Bellman- Gronwall lemma; we obtain 

\\v k ,2(t) ll<r-(exp(efc B t)-l)i (165) 

this shows that there exists 7 such that |u| 2 | < R for t < T e , which means 
that it is in L°°(0,T e ) for T e — ^; also, we have vj- in C(0,T e ) then as vj~ 
is solution of (|144l) . it is also bounded in C 2 (0,T £ ). 

□ 

Theorem 5.1. ( of Poincare-Lyapunov, for example see [12]) Consider the 
equation 

x — (A + B(t))x + g(t, x), x(t ) = x , t > t Q 

where x,Xo € R™, A is a constant matrix n x n with all its eigenvalues with 
negative real parts; B(t) is a matrix which is continuous with the property 
lim t ^ +oc ||£>(f;)|| = 0. The vector field is continuous with respect to t and x 
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is continuously differentiable with respect to x in a neighborhood of x = 0: 
moreover 

g(t,x) — o(||x||) when \\x\\ — > 
uniformly in t. Then, there exists constants C, to,#,/i such that if 

lkl|<C||xo||e-^- to ),t>to 

holds 
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